base_dir = "/Users/pol/owncloud/HiWi_folder/02_running_projects/02_MA_thesis/01_Data/Preprocessed/"
pt_path = paste0(base_dir, "PitchTiers/08_manual_corrected/")
Fujisaki_path = paste0(base_dir, "Fujisaki/")
library(contouR)

Linear regression slope compared

Compute slopes by drawing a linear regression line. Let’s first get an intuition by looking at the first 5 slopes.

files = list.files(pt_path)
files = files[grepl(".PitchTier", files)]
compute_regression_slope(files[1:5], pt_path, plot_only=TRUE)
## Loading required namespace: hqmisc

Now compute the slopes for all files:

features = compute_regression_slope(files, pt_path)
head(features)
##           name        a           b reg_RMSE naive_intercept  naive_slope
## 1   DF_ANG_VX1 26.14186  -9.7273180 2.572510        23.25355  -9.93577139
## 2 DF_ANG_VX10c 17.91798   0.1759461 3.313159        17.19068   2.13530953
## 3 DF_ANG_VX11a 22.49966  -3.4898652 3.705590        18.66072  -2.47602240
## 4  DF_ANG_VX12 21.70759  -2.1948611 2.689266        15.86251   0.05407044
## 5 DF_ANG_VX13a 28.59098 -10.0375849 3.931166        24.21989 -11.33068258
## 6  DF_ANG_VX14 20.33091  -4.5572098 4.454911        20.07695  -8.69039835
##   naive_RMSE min_max_slope min_max_RMSE
## 1   3.948228     -17.84009     4.188648
## 2   3.485802     -57.95575    28.774414
## 3   4.762231     -25.15022    14.607542
## 4   4.776461     -32.97042    18.349630
## 5   6.796655     -17.42301     5.193256
## 6   5.479758     -17.08930     8.409380

Let’s check the correlations. The regression intercept correlates strongly with the naive intercept. There is also a smaller positive correlation between naive and regression slope. However, it is not correlated with the slope between global minimum and maximum.

library(corrplot)
## corrplot 0.84 loaded
interesting_columns = 2:9
interesting_columns = names(features)[interesting_columns]
cor_df = features[, interesting_columns]
print(corrplot(cor(na.omit(cor_df)), method="number", type = "upper"))

##                            a           b    reg_RMSE naive_intercept
## a                1.000000000 -0.41068879  0.46425991      0.86460664
## b               -0.410688789  1.00000000 -0.03463125     -0.18848927
## reg_RMSE         0.464259908 -0.03463125  1.00000000      0.42803043
## naive_intercept  0.864606643 -0.18848927  0.42803043      1.00000000
## naive_slope     -0.487538388  0.53115991 -0.41912562     -0.54656440
## naive_RMSE       0.418484349  0.01567463  0.86820015      0.33379923
## min_max_slope    0.003059545  0.09329694 -0.06747011     -0.05259831
## min_max_RMSE     0.110679211  0.09607382  0.12344647      0.05688564
##                 naive_slope  naive_RMSE min_max_slope min_max_RMSE
## a               -0.48753839  0.41848435   0.003059545   0.11067921
## b                0.53115991  0.01567463   0.093296944   0.09607382
## reg_RMSE        -0.41912562  0.86820015  -0.067470114   0.12344647
## naive_intercept -0.54656440  0.33379923  -0.052598310   0.05688564
## naive_slope      1.00000000 -0.40522122   0.188722278   0.07816047
## naive_RMSE      -0.40522122  1.00000000  -0.031021720   0.13827723
## min_max_slope    0.18872228 -0.03102172   1.000000000   0.92176783
## min_max_RMSE     0.07816047  0.13827723   0.921767829   1.00000000
features$filename = features$name
features = add_meta_data(features)

interesting_columns = c("a", "naive_slope", "b", "naive_intercept")
for (col in interesting_columns){
  print(significance_test(na.omit(features), col))
}
## Loading required package: ggplot2
## Warning: package 'ggplot2' was built under R version 3.5.2
## Loading required package: magrittr

Phrase command

From these optimal hyper parameter combinations we can compute features on the phrase command

top_scores = read.csv("top_scores.csv")
phrase_features = compute_phrase_features(top_scores, Fujisaki_path)
names(phrase_features) = c("name",  "num_phrases", "Ap", "phrase_T0", "phrase_reg_slope", "phrase_reg_RMSE", "phrase_naive_slope", "phrase_naive_RMSE")
head(phrase_features)
##           name num_phrases      Ap phrase_T0 phrase_reg_slope
## 1   DF_ANG_VX1           1 0.73900     -0.46     -0.003836904
## 2 DF_ANG_VX10c           2 0.07395        NA               NA
## 3 DF_ANG_VX11a           1 0.25830     -0.47     -0.001377021
## 4  DF_ANG_VX12           1 0.28640     -0.46     -0.001097925
## 5 DF_ANG_VX13a           1 1.18280     -0.34     -0.005092311
## 6  DF_ANG_VX14           1 0.65570     -0.47     -0.003469269
##   phrase_reg_RMSE phrase_naive_slope phrase_naive_RMSE
## 1     0.009014474      -0.0034409529       0.022443360
## 2              NA                 NA                NA
## 3     0.034902102      -0.0015639914       0.079482718
## 4     0.005779488      -0.0009455769       0.009563485
## 5     0.038778261      -0.0040847909       0.067799577
## 6     0.022504489      -0.0031715768       0.038879622
features = combine_features(features, phrase_features, key = 'name')

We see that both slope computations (and also the RMSE) are highly correlated. Also we find a negative correlation between the amplitude and the slope

corrplot(cor(na.omit(phrase_features[, 3:8])), method="circle", type = "upper")

An visualize the features in boxplots:

phrase_features$filename = phrase_features$name
interesting_columns = names(phrase_features)[c(2:5, 7)]
for (col in interesting_columns){
  print(significance_test(na.omit(phrase_features), col))
}
## Warning: Computation failed in `stat_compare_means()`:
## data are essentially constant

ICCs

We can also devide the contour in it’s ICCs. Let’s first give some examples for the ICC separation:

compute_EAC_ICCs(files[1:5], pt_path, plot_only=TRUE)
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

We can now extract the features from the ICC as Mirjana Rajkovic, Jovicic, Grozdic, Zdravkovic and Subotic (2018) propose

feature = compute_EAC_ICCs(files, pt_path)
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in compute_EAC_ICCs(files, pt_path): No ICC detected
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in compute_EAC_ICCs(files, pt_path): No ICC detected
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped

## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
head(feature)
##           name num_ICCs num_pos_ICCs num_neg_ICCs mean_delta_F0
## 1   DF_ANG_VX1        3            0            3      6.763778
## 2 DF_ANG_VX10c        4            1            3      5.867933
## 3 DF_ANG_VX11a        6            1            5      6.336774
## 4  DF_ANG_VX12        6            1            5      6.107005
## 5 DF_ANG_VX13a        3            0            3      9.458832
## 6  DF_ANG_VX14        4            0            4      8.457135
##   sd_delta_F0   mean_tc      sd_tc mean_delta_tc sd_delta_tc         EAC
## 1    1.120784 0.2214286 0.14146645     0.2008929  0.12259159          NA
## 2    4.720117 0.2544643 0.20634854     0.2057143  0.20932654 0.003093019
## 3    4.255035 0.1767857 0.14047569     0.1545918  0.14104075 0.070351529
## 4    1.638377 0.1732143 0.08474903     0.1836735  0.08216504 0.053752638
## 5    7.077855 0.4535714 0.42077734     0.3535714  0.39753706          NA
## 6    6.395935 0.2303571 0.26455103     0.1885714  0.24742810          NA
features = combine_features(features, feature, key = 'name')

We can have a look at the correlations:

# Many
feature$filename = feature$name
cor_df = feature[, names(feature)[2:10]]
corrplot(cor(na.omit(cor_df)), method="circle", type = "upper")

cor_df = feature[, names(feature)[2:11]]
corrplot(cor(na.omit(cor_df)), method="circle", type = "upper")

An visualize the features in boxplots:

interesting_columns = names(feature)[2:11]
for (col in interesting_columns){
  print(significance_test(na.omit(feature), col))
}

Last but not least, we have to save the features in a CSV for later use:

write.csv(features, "slope_features.csv", row.names = FALSE)